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Abstract: Wc study quantum frequency translation and two-color 
photon interference enabled by the Bragg scattering four-wave mixing 
process in optical fiber. Using realistic model parameters, we compu- 
tationally and analytically determine the Green function and Schmidt 
modes for cases with various pump-pulse lengths. These cases can 
be categorized as either "non-discriminatory" or "discriminatory" in 
regards to their propensity to exhibit high-efficiency translation or 

C^ ' high-visibility two-photon interference for many different shapes of 

input wave packets or for only a few input wave packets, respectively. 
Also, for a particular case, the Schmidt mode set was found to be 
nearly equal to a Hermitc-Gaussian function set. The methods and 
results also apply with little modification to frequency conversion by 

^ sum-frequency conversion in optical crystals. 
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1. Introduction 

Quantum frequency translation (QFT) is the process whereby two spectral modes (or 
narrow bands of modes) are swapped in the sense that the properties of their quantum 
states (other than their central frequencies) are interchanged, as in Fig.[T]riJ[3J[3]. Such 
a process, of course, requires an external source or sink of energy, which is provided by 
a strong external optical field. Ideally, an arbitrary state of the narrow band of green 
modes will be swapped with the state of the narrow band of blue modes. These could be 
coherent, squeezed, or number states, and may also be entangled with other, separate 
degrees of freedom, such as other modes (e.g., violet and yellow) or states of atoms, 



quantum dots, etc. When the states of the green and blue modes are swapped, the prior 
entanglement should be preserved. This paper focuses on the QFT of single-photon 

states. 





Fig. 1. Quantum frequency translation of wavepacket states between green and 
blue spectral regions. The states being translated could be coherent, squeezed, or 
number states. 

QFT is defined to be "noiseless," that is, free of excess background noise. For example, 
if both modes to be swapped are initially vacuum, they should remain vacuum— there are 
no spontaneous processes that can populate modes in this case. The noiseless property 
distinguishes QFT from other frequency conversion processes such as optical parametric 
amplification, which has been widely studied for use in classical telecom systems [5J [3J 
[6]. In optical parametric amplification, power is spontaneously transferred from the 
pump(s) to the weak signal fields, and without any blue input power, the signal fields 
can build up spontaneously. 

The QFT capability will play important roles in quantum information science (QIS), 
where states of a quantum memory may be transmitted to a spatially remote quantum 
memory by emitting and absorbing photons. If the two memories operate using distinct 
frequencies of their absorption/emission lines, it is necessary to "translate" the commu- 
nication photons between these frequencies. Such a translation capability also allows 
sending the photons through low-loss optical fiber at wavelengths near 1500 nm, even 
in cases where the memories operate at other wavelengths. 

Two methods for all-optical QFT have been demonstrated: three-wave mixing in 
a second-order nonlinear optical medium (crystal) [2j [7] and more recently four-wave 
mixing (FWM) in a third-order nonlinear optical medium (silica fiber) [5]. In the case 
of three- wave mixing and referring to Fig. [TJ the strong pump field must have frequency 
(Op = (Ob — (O g . The configuration for FWM with two pumps, denoted p and q, is that of 
dynamic Bragg scattering, in which for example, one photon is removed from pump p 
and one is added to pump q. At the same time, one photon is removed from the green 
signal and one is added to the blue signal so that (O p + (O g = (O q + (Ob- This contrasts 
with optical parametric amplification, wherein photons are removed from the pump(s) 
and photons are added to both signals. 

Experiments have included QFT of one-half of a number-correlated state [7] and 
several examples of QFT of single-photon states by three- wave mixing in crystals [TJUJIl] 
and by FWM in fiber [8] . The main advantage of the FWM method over the three- wave 
mixing method is that it allows QFT between spectral bands separated by large or small 
shifts, whereas in practice three-wave mixing allows QFT only between far-separated 
bands [8]. This paper treats the theory of FWM, but the same methods can easily be 
applied to three-wave mixing. 

A potentially revolutionary application of QFT is that it is predicted to enable the 



interference of two photons (or other states of the modes) with different colors [TO] . 
This phenomenon generalizes the famous Hong-Ou-Mandel (HOM) interference effect, 
wherein two photons of like color and identical wave-packet structure impinge simul- 
taneously on a standard beam splitter. Because of the boson nature of photons, the 
beam-splitter output modes both contain either zero or two photons, but never one 
photon in each. This remarkable effect (which is analogous to the exchange interaction 
for electrons) enables a class of quantum information processing known as linear-optical 
quantum computing [11], including many variants [12] . 

The promise of QFT is that it offers the possibility to perform linear-optical quantum 
computing "over the rainbow." That is, one could use the light's color as a degree of 
freedom identifying optical qubits and perform linear-optical quantum computing using 
the generalization of the HOM effect involving distinct frequencies. This would entail 
creating entangled states involving many distinct frequencies (perhaps by methods such 
as in [13l E]), and manipulating such states using two- mode interference via QFT. 
This would require QFT between modes nearby in frequency, and so FWM would 
be greatly preferred over three-wave mixing. Other tools proposed for linear-optical 
quantum computing "over the rainbow" are temporal-mode filters ("pulse gates" [15] 
which could be used as channel add/drop filters in an optical communications system 
[El [17]) an d single-photon pulse shapers [HJE3], based on QFT by three- wave mixing. 

In this paper we extend our earlier treatments [3] [TU] to enable the explicit modeling 
of QFT by FWM in optical fiber or by three-wave mixing in a crystal. The treatment 
includes analytical and numerical components, and allows us to model realistically the 
effects of finite pulse durations, linear dispersion, finite phase-matching bandwidth, and 
complicating effects such as self-phase modulation (SPM) and cross-phase modulation 
(CPM) [20]. Previous theories of QFT have either assumed single- mode fields [3[, disal- 
lowing rigorous treatment of pulse propagation, or assumed that the signal pulse shapes 
are not altered significantly by the translation process [T5] [18] [19] , obscuring the accu- 
rate design of high-fidelity devices. 

The main results of this study are the following: We find that for realistic fiber 
and laser pulse parameters, essentially 100 percent efficient QFT can, in principle, be 
achieved. We show that for given temporal shapes of the pump pulses, usually there 
are unique shapes that the two signal wave packets should have in order that they 
be swapped optimally by QFT. These shapes are found as the optimal modes of a 
singular-value (Schmidt) decomposition of the Green function for the QFT process, 
and the shapes change significantly with changing translation efficiency. This result 
provides a method for designing QFT and two-photon interference with high efficiency 
and high fidelity. We also show, for a given pump temporal shape, how to design signal 
pulse shapes to create perfect HOM interference via the FWM process, as first predicted 
in [U3]. 

2. Theory 

The physical mechanism used to carry out frequency translation in fiber is the Bragg 
scattering (BS) four- wave mixing process J4][8][21], pictured in Fig. [2] For convenience, 
we refer to the two weak signal fields as "green" and "blue," although in principle their 
actual frequencies could correspond to visible or near-IR radiation. Two strong pump 
fields drive the process in which one photon from one of the pumps, pump p in Fig. [2] 
and one photon in the green field arc annihilated while one photon from the other pump 
(denoted q) and one photon in the blue field are created. The process conserves photon 
number in accordance with the Manley-Rowe relations [22| and can be summarized as 



Yp + Yg —> Yq + Yb, where y represents a photon. 
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Fig. 2. BS process where the pumps have frequencies (Op and CO q and the signals 
have frequencies co g and 0)/,, where 0)^ + (O g = (O q + CO/,. The arrows symbolize either 
annihilation (down) or creation (up) of a photon in that mode. The frequency (Oa 
is the average value of all the frequencies and Aco is the frequency separation 
between the pump fields and also between the signal fields. 

Before describing the realistic continuous-frequency case, we discuss a much simpler 
model in which each held is represented by a single- frequency mode. Such a model might 
apply to hclds in a cavity, but strictly speaking not to the pulse propagation scenario in 
which we are interested. The simple model has the advantage of permitting analytical 
solutions, thereby providing some insight into the QFT process [3j. For cases where 
each held can be described by a single-frequency mode, the Hamiltonian governing the 
process is given by 

,yf = A/3 (a g a g — a b a.b ) /2 + Ka g ab + K* a \a g , ( 1 ) 

where the creation/annihilation operators satisfy the usual commutation relations, and 
A/3 and K arc parameters quantifying the effects of phase- mismatch and nonlinearity [3] . 
The evolution is given by the spatial cquations-of-motion (with a minus sign compared 
to the usual time-domain Heiscnbcrg equation) 

-^ a b( s )(z)=i[a b ( g )(z),J^}. (2) 

The solution of ([2]) for a g iy) (z) as a function of position along the hber z for a fiber 
of length L is given by 

a g (z)=n(z)a g (0) + v(z)a b (0), (3) 

a*(z)=M*(z)a*(0)-v*(z)a,(0), (4) 

where the "transfer" functions are given by 

H{z) = cos(fe) + iAf5sm(kz)/2k, (5) 

v(z) =iKsin(kz)/k, (6) 

where k— (Aj3 2 /4+ Ik - ! 2 ) 1 / 2 . The transfer functions have the relation |^i| 2 + |v| 2 = 1, 
which ensures the process is unitary and conserves photon number. Mathematically, 
the BS translation process having frequency input and output "ports" is identical to the 
normal beam-splitter operation on degenerate monochromatic fields |23j . For perfect 
translation from green to blue, and vice versa, it must be that |v(L)| = 1 for a medium 
of length L. In this case the quantum states of the modes are completely swapped. For 
example, a single green photon is replaced by a single blue photon, that is, 1 1 , 0) — >• |0, 1) , 



where \i,j) denotes the state in which i photons are in the green region and j photons 
are in the blue region. 

Since the BS process is analogous to beamsplitter action, all of the effects that a 
beamsplitter can induce on a field will also occur in the BS process. In particular, the 
effect analogous to two-photon HOM interference will occur, where two independent 
photons of differing frequencies will interfere with one another and give a perfect HOM 
dip under the appropriate circumstances [10| . Just as in the monochromatic HOM inter- 
ference effect, the "frequency beamsplitter" must have equal transmission and reflection 
coefficients, that is \jx\ 2 = \v\ 2 = 1/2, so that it is equally likely for the blue (green) pho- 
ton to stay in the blue (green) mode or be translated to the green (blue) mode. Using 
the above equations it can be shown that in this case, for an input state consisting of 
one green photon and one blue photon, denoted by a^(0)a' b (0)\vac) = |1, l)m, the output 
state is given by 

|l,l) m =f[|Ai| 2 -|v| 2 ]|l,l> out +\^Aiv|2,0) out -V2^v*|0,2) out N ). (7) 



For equal values of \jJ,\ 2 and |v| 2 the |1, 1) component of the state cancels to zero. 

The above analysis is valid only if all the fields involved are single- mode. For fields that 
contain multiple-frequency components, as any pulsed field does, the theory becomes 
much more complex. For arbitrary input, and even "simple" multi-mode input such as 
Gaussian pulses, the output cannot in general be determined analytically. However, it 
is possible to assign analytically a Green function G in either frequency (0 or time t that 
relates the input wavepacket to the output wavepackct for arbitrary input, analogous 
to ([3]) and ([3]), for given pump spectra, fiber dispersion, and fiber length [10]. Denoting 
the creation operator at position z by a'(z,(o), the input a'(Q,(o), and output a'(L, (o) 
are related by the forward transformation 



a f (L,0))= d(0G* (0,(8^(0, to). (8) 

For our purposes, the backward transformation is more useful: 

a f (0,co)= /<ico'G(a),a)')fl t (L,a)'), (9) 



which gives the input operator in terms of an integral over the output operator. 

A single-photon input wavepacket given by J d(OA(0, <3))a^ (O,co)\vac) , where A(0,<») is 
the spectrum of the input wavepacket, is translated in the manner given by ©. The 
state can thus be expressed in terms of the output operator as 



\W)= dcoA(0,co) / deo G((o,G> )a r (L,eo )\vac). (10) 

The Green function is unitary and therefore satisfies the condition 

d(oG((o,(o')G*((o'' ,ca) = 8((o-(o"). (11) 



To facilitate an intuitive understanding of the transfer function, we break up frequency 
space into two pieces; the range that contains the green and the range that contains the 
blue. Each operator or variable will be denoted by a "g" or "£>" subscript. For example, 
the creation operator for an green photon at a particular green frequency will be written 
as a'Jco s ). In this case the transfer function can be written in block form as |10j 



G gg {a g ,G) g ) G gb ((O g ,(O b ) 
G bg (a> b ,a> g ) G bb (co b ,co b ) 



(12) 



where G xy refers to the evolution of an output creation operator in range y which came 
from range x. 

Although the G function completely describes the evolution of arbitrary inputs, it 
is difficult to discover which input pulses give rise to desirable effects such as high- 
efficiency translation or good two-color HOM interference just by examining the G 
function. However, by performing a Schmidt (singular-value) decomposition (SVD) [24] 
on the G function or on the sub- matrices G. rv , the input wavefunctions necessary to 
achieve these effects can be found [10]. For example, consider the translation of an 
input green field to the blue field, which is described by the Green sub- matrix by G g b, 
which is not Hcrmitian and therefore does not admit a unitary decomposition. We write 
Ggb in its SVD form as 

G gb (co,co') = -£p„y„(fi))w*(fi)'), (13) 

where the p„ are the real, positive singular values (the generalization of eigenvalues), 
and V„ and w„ are the «th column vectors of the unitary "matrices" V and w describing 
the decomposition, respectively. The most useful way to interpret this decomposition 
is that a given green wavepacket V„ will be translated to the blue wavepacket w„ with 
probability p^. Hence, only input green wavepackets V„ that have associated singular 
values p„ = 1 will translate over to the blue mode with 100 percent probability. The G gg 
"matrix" describes an input green wavepacket scattering within the green spectral region 
at the output (no translation from green to blue, although the shape of the wavepacket 
may change). If the SVD of G gg is given by [TU] 

G gg {(0,(0)=^ n V n {(0)V* n {m), (14) 

n 

then T„ is the probability that an input green wavepacket V n will scatter to an output 
green wavepacket v„ (the V matrix here being the same as in the G g b decomposition). 
The singular values satisfy the relation rf + p„ = 1 , as they must to conserve photon 
number. 

Similarly, the Gb g and Gbb have their own SVDs, given by 

G bg {co,(o')=Y,PnW„((o)v*(co') 1 (15) 

n 

G bb {(o,(o)=Y i 'C n W n {(o)w* n {(o), (16) 

// 

where the w and v matrices, and the singular values T and p, are the same as in (fl~3|) and 
([j"4]) . This highlights the interconnectedncss between the various sub-Green functions 
in that there are only four unique Schmidt modes for each index n that completely 
describe the process. Before treating exclusively single-photon inputs, we emphasize 
that the Green functions determined by our present approach can be used to treat 
arbitrary multi-photon inputs as well. 

The SVD technique can be used to discover what input wavepackets are suitable for 
optimal two-photon HOM interference |10j . The condition is somewhat intuitive. Since 
in the single-mode, two-color case it is necessary for both green and blue inputs to have 
a 50 percent probability of translating to the other mode for perfect HOM interference, 
it is reasonable to guess this condition is also true in the multi-mode two-color case. 
Hence for the cases in which x^ = p% = 1/2, the corresponding input vectors V„ and 



W„ will show perfect two-color HOM interference. Alternatively, this condition can be 
written as 

0-„=2T„p„ = l, (17) 

We refer to the quantity o~„ = 2x n p n as the HOM singular value, because it appears 
as such for a Schmidt decomposition kernel designed specifically to optimize the HOM 
interference. (See Eq. (37) of [TO]-) The closer it is to unity, the better the HOM inter- 
ference. 

There is an alternative, complimentary, method to calculate the degree of HOM 
interference that a BS process will exhibit, for which the Green function does not have 
to be calculated. First, consider the case in which there are two photons, either of which 
could be in the green or blue region, or in some superposition of both regions. In general, 
the state of this bi-photon can be written as 

|V/)=C 11 |l,l)+C 20 |2,0)+Co 2 |0,2). (18) 

The entire state |y/) must be normalized and the component states are orthonormal. 
Now, for input consisting of one green photon and one blue photon having input spectra 
of A g (0,(Og) and At, (0, (Ob) respectively, the state is 

dcO g dC0bAg(O, O) g )A b {0, (Q h )al(0, co ? )a£(0, (Q b )\ vac) . (19) 

If the photons undergo the BS process, the input photon operators evolve according to 
©. By switching the limits of integration and noting that [J dcoAx(0,Co)GxY{(O,(>h')] - > 
Axy{L,<X>y), where both X and Y could be either g or b and Axy{L,(Oy) is a component 
of the spectrum at the output of the fiber, the state can be written compactly as 

\W)=[ dcOga\(L, (Og)Agg(L, (O g ) + / d(O s a\{L, (o b )A bg {L, a> b )] 

r , , r , , ( 2 °) 

x [ / dco' g al(L, co'g)A gb (L, aj g ) + / dco' b al(L, co' b )A bb (L, co' b )] \ vac) . 

The probability, given a particular output state, for the various sub-space states to 
occur is given by the modulus squared of the coefficients of these states, |C, ; | 2 . For the 
|1,1) state the coefficient is found by acting the aforementioned operators, with their 
accompanying functions, on the vacuum, yielding 

IVn) = / / d(Ogd(O b [Agg(L,co g )A bb (L,co b ) +A gb (L,cOg)A bg {L,co b )] \lm g ,lco b ) (21) 

The first term, Ag g (L, (O g )A bb (L, co b ), corresponds to the case where a photon initially 
in the green mode (i.e. wavepacket) scatters within the green mode at the output, and 
a photon initially in the blue mode scatters within the blue mode at the output. The 
second term, A gb (L, 03 g )A bg {L, (O b ), corresponds to the case where a photon initially in 
the green mode scatters to the blue mode at the output, and a photon initially in the 
blue mode scatters to the green mode at the output. Since these terms are being added 
before their modulus squared is taken, there is the chance that they could cancel if the 
terms are proportional and out of phase with one another. This is exactly what happens 
in the case in which ideal HOM interference occurs, leaving zero probability amplitude 
for the field to emerge in the |1,1) state. The expression for the probability for this 
state to occur P\ \ is 

Pn= J J dco g d(O b [\Agg{L,cOg)\ 2 \A bb {L,co b )\ 2 +\Ag b {L,co g )\ 2 \A b g{L,co b )\ 2 + 

2Re[A g g{L, (Og)A bb {L, (O b )A* gb {L, (Og)A* bg (L, co b )\}] . 



This gives an expression to calculate the coincidence probability for any given input 
pulses without knowing the Green function for the process. In the alternative case in 
which we launch particular Schmidt modes as input signals, we showed previously that 
the coincidence count rate equals P\\ = (tf — p 2 ) 2 [TO]. Therefore, in this case we have a 
simple relation between fn and the HOM singular values; Pu + a 2 = 1. Other important 
quantities are P20 and P02, the probabilities that two photons appear in one or the other 
output ports together. The expressions for these, and their derivation can be found in 
Appendix A. 

3. Numerical implementation 

Unfortunately, it is not possible to discover the Green function and its Schmidt modes 
by examining the Hamiltonian directly. One must solve for the full propagation first, 
then examine the input-output relations. This emphasizes that the Schmidt modes are 
global properties, depending on the entire evolution. For example, the input Schmidt 
modes depend on the length of the medium. 

The appropriate Hamiltonian generates the well-known nonlinear Schrodingcr equa- 
tion, which needs to be solved. In most realistic circumstances, taking into account 
actual field shapes and real fiber dispersion properties, it is not possible to solve the 
nonlinear Schrodingcr equation for the fields analytically. So we turn to numerical com- 
putation to determine the evolution of a completely specified input state to a completely 
determined output state. The key difficulty with such an approach is that, in general, 
nonlinear operator equations cannot be simulated by any simple numerical scheme. 
When the evolution is both in space and time as in our case, the problem becomes one 
suited for techniques such as the positive-P representation [25l [26] . Fortunately, in our 
case, although the evolution is nonlinear in the pump fields, it is linear in the weak sig- 
nal fields. This allows a simplification if we treat the pump fields as classically evolving 
parameters in the quantum Hciscnbcrg equations of motion for the weak signals. Then 
the field operators for the weak signals can be solved in terms of Green functions that 
obey the classical evolution equations. Therefore, we solve numerically for the classi- 
cal evolution of the Green functions and use those solutions to calculate any quantum 
probability or correlation function concerning the signal fields. An analogous procedure 
was used in [2"71l2"5] . 

The usual starting point for solving for the nonlinear evolution of a pulse is to write the 
total field as E(z,t) = A(z,t)e'^P° z ~ C00 '' , where A(z,t) is a slowly varying envelope function, 
0>o is the high frequency of the carrier wave, and j3o is the propagation constant at do- 
In that form the evolution is governed by the nonlinear Schrodinger equation 

d z A{z,t) = iP(id t )A(z,t) + iY\A(z,t)\ 2 A{z,t). (23) 

Although this is effective for single pulses, for electric fields made up of multiple, distinct 
pulses with significantly different frequency components, solving this form of the equa- 
tion is unnecessarily computationally intensive and impractical. Instead, we solve four 
coupled equations, derived from the nonlinear Schrodinger equation, in which terms re- 
lating to dispersion, self-phase modulation, cross-phase, and the BS process were kept. 
Numerically, the frequency components of each field are associated with a frequency 
mesh pertaining only to that field. This drastically reduces the number of calculations 
needed to compute the solution in that it is unnecessary to resolve the large empty band- 
width in between fields. In this scheme, the total electric field is the sum of the four 
individual fields, each having a slowly varying amplitude function Aj centered around 



a carrier wave of frequency COj; 

E(z,t)=Y j A J (z 1 t)e i ^-^+c.c. : (24) 

j 

where j denotes p, q, g or b. In quantum theory, A ; - is proportional to an annihilation 
operator. Substituting this ansatz into the nonlinear Schrodinger equation leads to, 
with suppressed independent variables z and t, the four coupled BS equations (where 
d z represents the partial derivative with respect to z) 

d-Ap = ij3 p {id T ;co p )A p + I /y|A ; ,| 2 + 2/y£ \A k \ 2 \A p + 2iyAlA q A b exp[-iAI5z} 



(25) 



d z A q = ij5 p (id T ;co q )A q + jy|A ? | 2 + 2r'y£ |A A .| 2 A <? + 2/yA^4 p A^exp[/A/3z] 

V **? J 

d-A x = ip p (id T ;C0g)A g + ( /y|A/ + 2/y£ |A*| 2 j A ? + 2/yA; ;i A c/ A i exp[-/A/3z] 

d-A h = ip p (id T ;cQ b )A b + ( /y|Ai| 2 + 2/y£ |A A .| 2 ) A 6 +2ryA^A / A g exp[rAj3z], 

V W* ) 

where A/3 is the dispersive mismatch given by Aj8 = j3 ; , + J3 S — (f5 q + j3^) + y(P q — P p ), P q 
and P p are the peak pump powers, and y is the nonlinear coefficient. The first term 
on the right hand side describes the effect of linear dispersion, the next two terms 
describe SPM and CPM respectively, and the last term describes the BS four-wave 
mixing process. The quantity j3 p (idT\£lo) is given by 



w » ;ao ),£ ggslW! _p,.) W ,.». 



(26) 



n=\ 



where )3'"'(i2o) is the nth derivative of /3 evaluated at Q.q, and T = t — j5^ l '(a)p)z- In this 
scheme all the fields are propagating in the arbitrary choice of the frame of reference of 
pump p. With one frame of reference there need only be one universal time mesh. 

Although (|2"5|) is for quantum field operators, we can realistically replace the pump 
fields A p and A q by classical functions if these fields are in strong coherent states. We 
solve the first two equations for the pumps ignoring the small effects of the weak signal 
fields A g and A&, but retaining the SPM and CPM effects of the strong pumps. We 
then use these pump solutions in the second pair of equations for the weak signal fields. 
Because these equations are linear in the operators, we can treat them formally as 
classical equations. Then the quantum effects can be fully accounted for by using the 
commutation relations when calculating quantities such as Pn, as describe above. 

To solve the coupled equations numerically we use the split-step Fourier method 
(SSFM) [29l |30]. For each electric field Ej we approximately solve the operator equa- 
tion d z Aj(T,z) = (Dj +Nj)Aj(T,z), where Dj and Nj are the dispersive and nonlinear 
differential operators, acting in the frequency and time domains respectively, which act 
on Ej given in Eq. (f2~5j) . The approximation is that A(T,z) is taken to A(T,z + Az) by 

A(T,z + Az) « rap (—£ leap (/"* ' Ndzjexp(^DjA(T,z). (27) 



where Az = L/N and N is the number of steps. The nonlinear component integral is 
found by use of a fourth-order fixed-step Runge-Kutta method. 

The SSFM allows one to calculate a particular output amplitude for a particular 
arbitrary input amplitude. However, discovering the Green function that governs the 
process requires solving the propagation equation for a family of orthogonal input am- 
plitudes. Let A,„ be an n x n matrix with each column equal to one time-domain vector 
from an orthogonal set, A ou , be an n x n matrix with each column equal to the out- 
put of the n column of A/„, and G be the Green function. In this paper we use the 
Hermite- Gaussian (HG) functions as the orthogonal set, where the n th function is given 
by 

Wn(x)=( ^) e-- v2 / 2 (-l)V v2 — e- x \ (28) 

Y V ' \2»n!^/n) V ' dx" v ' 

where the characteristic time of the basis set is set by the overall time window of the 
mesh. Due to the linear nature of the Green function formalism, as expressed in ©, 
the forward Green function expressed in matrix notation is 

G^=A out -Ar n \ (29) 

This is strictly true only if all of the output vectors of A out can be described as some 
linear combination of the orthogonal vectors of A,„. In practice, an orthogonal vector 
of order j will evolve to have components of both lower and higher order than j. In 
the cases where the yth mode evolves to include components of higher order than «, G 
won't satisfy the unitary conditions of (TTTJ) , which is often the case for the evolution 
of the very-high-order modes. In that case, only the subset of orthogonal input vectors 
that evolve into output that can be described by all n input vectors can be considered 
"valid" output, as far as determining the Green function is concerned. Let A,' nsu £ be 
a m x m dimensional matrix that is the sub-matrix of A,„ where all the input vectors 
evolve to output vectors (described by the n x m matrix A outsu b, where m^n) that can 
be described as a linear combination of the vectors of A;„. Then the Green function that 
always satisfies the unitary conditions is given by 

G sub = A out,sub ■ A i, usu bi ( 30 ) 

and will be a matrix of dimensions n x m. In practice, the simplest way in which to 
determine the correct matrix subsets is by applying the unitary conditions (calculating 
G-G^) and seeing for what subset of vectors they hold. When that is determined, G* b 
can then be calculated. 

Finally, there is the question of fiber modeling. As these calculations are meant to help 
explain past experimental results [8] and guide future experiments using solid-core pho- 
tonic crystal fiber [31] , we desire a fiber model that fairly accurately describes photonic 
crystal fibers. We employ the often used step- index model in describing photonic crystal 
fibers [32j , which approximates the photonic crystal fiber as a step-index fiber having a 
core index n core (usually the index of bulk fused silica) and an effective cladding index 
n dad, where n c i ac /(a>) = /+ (1 — f)n C ore{(o). Here / is the air-filling fraction of the lattice 
region compromising the cladding, which along with the core radius a, parameterizes 
the fiber. By measuring the modulation instability wavelength signals produced by a 
fiber as a function of pump wavelength, best-fit values for / and a can be determined. 
For a particular fiber that has been used in past experiments and is a good candidate for 
future experiments, this procedure has been carried out, resulting in parameter values 
of / = 0.494 and a = 0.72 fim, with the resulting dispersion parameter D shown in Fig. 
|3l This is the dispersion profile used for all numerical calculations in this paper. 
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Fig. 3. The D parameter for the photonic crystal fiber used in simulations from 
600 nm to 1600 nm. The dashed line is at D = for reference. 



Given a particular dispersion profile of a fiber and the central wavelengths at which 
translation will occur, the fiber can effectively translate input signals only with spectral 
widths smaller than the characteristic bandwidth of the BS process. As seen in equations 
(O in the case of low conversion efficiency, the state is proportional to sinc(Aj3L/2), 
which defines the dispersive mismatch of the process. Figure [4] shows this sine function 
for the fiber, central wavelengths, and length (20 m) used in the simulations. All fields 
were taken to be monochromatic, and the pumps where fixed at 808 nm and 845 nm. 
The "blue" signal field was varied from its central wavelength of approximately 649 nm 
(shown in frequency on the horizontal axis) while the "green" field was varied from its 



central wavelength of 673 nm to conserve the energy of the process (co s 



(On 



(Oi-Oq). 



The full-width at half-maximum (FWHM) of the central sine lobe is approximately 
0.3 Trad/s. This number gives a rough estimate of the effective bandwidth of the BS 
process. As we will see, all else being equal, input signals that have spectral widths 
larger than this number do not translate as efficiently as signals that have spectral 
widths lower than this number. 
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Fig. 4. Phase-matching function sinc(A/3L/2) for the fiber used in the simulation, 
where the pump wavelengths were fixed and the "green" and "blue" fields fre- 
quencies were varied in an energy-conserving manner. The horizontal axis is the 
"green" frequency variation from its corresponding central wavelength value of 
673 nm. For this plot the fiber length is L = 20m, and the FWHM of the function 
is approximately 0.3 Trad/s (10 12 rad/s). 



4. Numerical results 

This section details the results of simulations for several different cases, broken up into 
two main categories; cases that highlight properties relating to effective total translation 
from one signal field to another, and cases that highlight properties relating to good 
two-color HOM interference. For the cases in which the Green function of the process 
was calculated, the two pumps were identical transform-limited Gaussian-shaped fields. 
The pump intensity profiles had either a "long" FWHM duration of 1000 ps, making 
them quasi-CW for most input signal fields used, or a "short" FWHM duration of 70 ps. 
At the input to the fiber, the centers of the fields were placed at zero time. The central 
wavelengths for the pump fields p and q were 808 nm and 845 nm, respectively, placing 
them in the anomalous-dispersion region. The signals "green" g and "blue" b were at 
wavelengths 673 nm and 649 nm, respectively, placing them in the normal-dispersion 
region. The 845-nm and 649-nm fields had nearly equal group velocities (= 2.0044 x 10 8 
m/s), and the 808- nm and 673- nm fields had nearly equal group velocities (= 2.0049 x 
10 8 m/s). The nonlinear coefficient 7 was set to 100 W^'km - for all calculations. 

The Green functions were calculated in the HG mode basis, where the input and 
output field amplitudes were decomposed into a common set of HG functions. The 
Green function for a particular case was found by sequentially setting the green or 
blue input field to a single Hermite-Gaussian function, solving the propagation, and 
projecting the output pulse onto the HG basis functions. Thereby the green function is 
represented in the HG mode basis. The characteristic time scale for the HG functions, 
which we define as the FWHM of the zeroth order HG function, used to produce the 
figures was found by first arbitrarily choosing a characteristic time of 40 ps, fitting 
the resulting absolute value of the first Schmidt mode to a Gaussian functional, and 
calculating the Green function again using this new characteristic time. This time-scale 
choice for defining an "optimal HG basis" produced Green function matrices that were 
the least complex, some of which were nearly diagonal, and therefore the simplest to 
understand. The HG basis set was centered in time at the peak of the pump and signal 
inputs. 

Parameter studies were also carried out in which the two pump pulses' intensity dura- 
tions were set at 1000 and 70 ps, and the input signal field durations were varied. Either 
the translation efficiency or HOM interference properties were documented, depending 
on the main category of the case. 

There are a few ways in which to gauge the accuracy of the numerical solutions. While 
not explicitly built into the calculations, for all calculations the sum of the green and 
blue pulse energies (calculated in the frequency domain) were found to be conserved, 
in accordance with the Manley-Rowe relation [22] for the BS process. Also, there are 
unitary conditions on the Green functions that govern the process, given by (|TT|) . Al- 
though checking these conditions is significantly more expensive computationally than 
verifying the Manley-Rowe relations, and therefore is not feasible for studies requiring 
large numbers of cases like the below- mentioned parameter studies, this check was car- 
ried out for several representative cases for such studies. The Green functions were then 
checked against (|lll) . and in these cases were found to agree well. Finally, the appro- 
priate number of steps were taken for each calculation so that the solutions converged 
and were not meaningfully dependent on the exact number of steps. 



1 In most cases the first Schmidt mode was quite Gaussian-like. 



4--1- Frequency translation, long pulse 

Figures [SJ [HI and [7] detail the case where the pump pulses were 1000 ps long and had 
peak powers of 400 mW. It was found by trial and error that translation occurred 
efficiently at a fiber length of 20 m for pumps with peak powers around 400 mW. It 
was found that the characteristic time of the optimal HG basis was w 243 ps. Figure 
[5ja) shows the absolute value of the V matrix, as defined in (|13|) , from an SVD of the 
G g b Green function. This Green function corresponds to the case of translation from 
the green mode to the blue mode. As discussed previously, the columns of the V matrix 
represent the input Schmidt modes, so in this case the V matrix represents the input 
signal field. That is, 

V„(co) = £c" jVj (co). (31) 

./' 

The horizontal axis is the Schmidt mode number n, where order was determined by 
the magnitude of the corresponding singular value p, largest to smallest, for the first 
25 Schmidt modes. The vertical axis is the HG mode number j of the Schmidt modes 
using the same HG basis used to find the Green function. The plotted quantity is \c"A. 
For example, the first Schmidt mode of the V matrix is mostly composed of the first 
HG function and almost completely described by the first four HG functions. Thus 
we can conclude that the temporal duration of the first Schmidt mode is close to the 
characteristic time of the optimal HG basis, 243 ps. 

The output Schmidt modes, which physically are the results of propagating each cor- 
responding input Schmidt mode through the QFT process in the fiber, are described 
by the w matrix. The absolute value of this matrix, the column vectors of which cor- 
respond to output in the blue mode, is shown in Fig. [5jb). The axes are the same as 
those in part (a). Interestingly, the V and w matrices are qualitatively similar, although 
the w matrix Schmidt modes tend to be dominated by either even or odd HG functions 
whereas the V matrix Schmidt modes are composed of a more or less equal number of 
even and odd HG functions. 

The black circles of the red line in Fig. [5jc) are the squared singular values p of 
the Ggb Green function for this case. The circles of the blue line correspond to the 
squared singular values T 2 of the G gg Green function. The p 2 describe the translation 
efficiencies of the input green Schmidt modes to the output blue Schmidt modes, whereas 
the T 2 describe the "non-translation" efficiencies of the input green Schmidt modes to 
the output green Schmidt modes . This corresponds to the transmission and possible 
reshaping of the green mode. For the Manlcy-Rowc relations to be satisfied, the square 
of these singular values must add to one, which they do numerically in this case (as well 
as in all other cases) to one part in one thousando As is evident from the figure, many 
of the modes translate efficiently with the first several having efficiencies of over 90 
percent. We call this configuration "non-discriminatory" in the sense that it effectively 
translates many different dissimilar modes efficiently. 

Figure [5] shows the translation efficiency as a function of length long the fiber for 
the first three green input Schmidt modes. Here translation efficiency is defined as 
the number of photons in the blue mode divided by the initial number of photons 
in the green mode. Note that there were initially no photons in the blue mode. The 
red, green, and blue lines correspond to the first, second, and third Schmidt modes, 
respectively. The efficiency for these modes is very sinusoidal, almost as if the signals 
were monochromatic and described by §5§ and ^. 



2 The higher-order singular values are accurate to one part in one thousand, but the lower-order 
singular values are usually more accurate, on the order of one part in one million. 
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Fig. 5. Singular values, V matrix (columns of which are coefficients in the HG basis 
for individual Schmidt modes), and w matrix of G„t, Green function relating to the 
case with long pumps (1000 ps) having peak powers of 400 mW and a fiber length 
of 20 m. (a) Absolute value of the V matrix (input Schmidt modes), (b) Absolute 
value of the w matrix (output Schmidt modes), (c) Squared singular values p 2 
(red) and T 2 (blue) indexed by the corresponding Schmidt mode number. 



Figure [7] shows the absolute amplitudes and phases in both time and frequency of the 
first three input Schmidt modes (relating to the V matrix) of the green (673 nm) mode. 
The input amplitudes of the pump fields are also shown as a black dashed line. As could 
be guessed from Fig.[5ja), the modes look very much like HG functions, although there 
is some asymmetry and delay in both time and frequency. These effects are likely the 
results of pump evolution, which includes convection, dispersion and CPM. Figure [TJa) 
shows that first few Schmidt modes are much narrower in time than the pump fields. 
Looking at Fig. EJb), which shows the phases in time, all the input Schmidt modes, 
disregarding the n phase jumps, are distinctly parabolic with positive curvature. For 
the phase convention we use, this shape implies the Schmidt modes at the fiber's input 
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Fig. 6. Translation efficiency as a function of length along the fiber relating to 
the case with long pumps (1000 ps) having peak powers of 400 mW. The red, 
green, and blue lines relate to the first, second, and third green Schmidt modes, 
respectively. 



are down-chirped (frequency decreases in time), meaning that when they propagate 
through the medium, which for the signals is normally-dispersive, they will temporally 
compress. That is, they are dispersion pre-compensated. CPM also chirps the signals 
in the same way as linear dispersion. Figure [7fc) shows that the pump fields are nearly 
monochromatic as compared to the Schmidt modes, although all modes are well within 
the translation phase-matching bandwidth of the fiber. 

4-2. Frequency translation, short pulse 

Figures 13 G3 and [TO] detail the case where the pump pulses were 70 ps long and had 
peak powers of 400 mW. The goal in studying this case in relation to the former case is 
to highlight how using temporally shorter and spectrally broader pump pulses impacts 
the translation properties of the process. It was found that the characteristic time of 
the optimal HG basis for this case was « 43 ps. Figure Eta) shows the absolute value 
of the V matrix (input green Schmidt modes) from the SVD of the G g y Green function, 
while Fig. E£b) shows the absolute value of the w matrix (output Schmidt modes) . The 
horizontal and vertical axes have the same meaning as they did in the previous case. 

While the absolute values of the V and w matrices are qualitatively similar, they 
are significantly different from their counterparts of the previous case. While they start 
out for low mode order approximately diagonal, like in the previous case, the Schmidt 
modes quickly become seemingly random and unordered. This is likely due to these 
Schmidt modes' singular values [the squared value of the singular values are shown on 
the red line of Fig. Etc)], which quickly drop to almost zero, signifying that almost no 
translation occurred. Since there is an infinite number of Schmidt modes with singular 
values near zero (practically any random combination of HG modes generally leads to 
a low singular value), even Schmidt modes of relatively low Schmidt mode number that 
have singular values close to zero can look random (non-uniquely defined) in the HG 
basis. Since in this case only the first few modes translate with significant efficiency, 
we call this configuration "discriminatory" in the sense that only a few closely related 
modes effectively translate. 

Figure |9] shows the translation efficiency as a function of length along the fiber for 
the first three input green Schmidt modes in the short-pumps case. The red, green, and 
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Fig. 7. Amplitude and phase in both frequency and time of the green 673-nm 
input Schmidt modes relating to the case with long pumps (1000 ps) having peak 
powers of 400 mW. The red, green, and blue lines relate to the first, second, and 
third Schmidt modes, respectively. Note: Trad/s indicates 10 rad/s. 



blue lines correspond to the first, second, and third Schmidt modes, respectively. Unlike 
the former long-pumps case, the efficiency curves of the modes are drastically different, 
relating to the fact that their singular values are quite different. 

Figure [10] shows the absolute amplitudes and phase in both time and frequency of 
the first three green 673 nm input Schmidt modes (relating to the V matrix). The input 
amplitudes of the pump fields are also shown as a black dashed line. Qualitatively, the 
amplitudes for this case are similar to the amplitudes in the former case; although asym- 
metric, they are roughly HG-likc. Quantitatively, these amplitudes are much broader 
spectrally than in the former case. The phases in time are essentially parabolas with 
positive curvature much larger than that of the previous case, as expected because the 
bandwidths are much larger. 

Figure [11] shows a parameter study where the pump pulses are identical Gaussians 
with temporal FWHM durations of 1000 ps and peak powers of 400 mW. In Fig. QT] 
(a) , the maximal translation efficiency achieved along the fiber is plotted as a function 
of input signal pulse FWHM duration. For relatively "long" pulses in low hundreds of 
picoseconds, the translation efficiency is quite high due to the near-CW nature of the 
pump and the width of the translation bandwidth. As the pulses become shorter, par- 
ticularly around 43 ps, which corresponds to a spectral bandwidth of 92 Grad/s, the 
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Fig. 8. Singular values, V matrix, and w matrix of Gob Green function relating to 
the case with short pumps (70 ps) having peak powers of 400 mW. (a) Absolute 
value of the V matrix, (b) Absolute value of the w matrix, (c) Squared singular- 
values p 2 (red) and T 2 (blue) index by the corresponding Schmidt mode number. 



pulse spectral width begins to match that of the translation phase-matching bandwidth. 
Therefore, this dramatic drop in efficiency is likely caused by the signal spectral band- 
width being too large for the BS process to accommodate. In Fig.QjJb), the length at 
which the maximal translation occurs is plotted as a function of input signal duration. 
The curve has the same qualitative shape as the curve in part (a); for short signal du- 
rations, the length at which maximal translation occurs is shorter than that of signals 
with long durations. This likely is due to the fact that short pulses have large spec- 
tral bandwidths, which because of CPM from the pumps, become even more spectrally 
broad as they evolve, making them in harder to translate. Hence the best translation 
occurs for these signals when they haven't evolved much. Signals with long duration 
and small spectral width do not spread out in frequency as much as the short signals 
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Fig. 9. Translation efficiency as a function of length along the fiber relating to 
the case with short pumps (70 ps) having peak powers of 400 mW. The red, 
green, and blue lines relate to the first, second, and third green Schmidt modes, 
respectively. 



do, so they are less limited by the magnitude of the translation bandwidth and achieve 
good translation along the entire fiber. 

Figure [T^] show a parameter study where the pump pulses are identical Gaussians 
with temporal FWHM durations of 70 ps and peak powers of 400 mW. The efficiency, 
shown in part (a), for short signal inputs behaves similarly to the previous long-pumps 
cases; around 20 ps the efficiency drops sharply due to the signal having larger transla- 
tion phase-matching bandwidth than the fiber. Unlike the previous case, the efficiency 
reaches a maximum around 35 ps and slowly declines for wider input signal widths. 
Relatedly, it is at about this signal width that the length at which the maximal effi- 
ciency occurs, shown in part (b), is equal to the length of the fiber, suggesting that 
the translation efficiency would be higher for a longer fiber or higher fiber for these 
signal widths. As the signal width becomes equal and then greater to the pumps widths 
there are two main elements affecting the efficiency; the lack of power/fiber length and 
the fact that the pump pulses no longer fully encompass the signal pulse in time. Both 
of these lead to lower translation efficiency. But it is important to note that for some 
signal width range, from about 25 ps to 60 ps, reasonably good translation efficiency 
will occur for short pump pulses with reasonable parameters such as peak pump powers 
of 400 mW and a fiber length of 20 m. 

4-3. HOM interference, long pulse 

Figures HH [HI and [15] detail the case where the pump pulses were 1000 ps long and 
had peak powers of 200 mW. This configuration was designed to exhibit good two-color 
HOM interference by lowering the pumps' peak power so as to lessen the translation 
efficiency, since the HOM interference effect is optimal when 50 percent, not 100 percent, 
translation occurs. It was found that the characteristic time of the optimal HG basis 
was w 146 ps. Figure flBTa) shows the absolute value of the V matrix, as defined in (fl"3"|) . 
from an SVD of the G g b Green function. This Green function corresponds to the case 
of translation from the green mode to the blue mode. Figure [TBTb) shows the absolute 
value of the w matrix. The horizontal and vertical axes are the same as they were the 
previous cases. Surprisingly, both the V and w matrices are nearly diagonal, meaning 
that the input Schmidt mode set and the HG basis functions are nearly identical. In 
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Fig. 10. Amplitude and phase in both frequency and time of the green 673-nm 
input Schmidt modes relating to the case with short pumps (70 ps) having peak 
powers of 400 mW. 



Appendix B wc present a theoretical model explaining this result. The long pump pulse 
of 1000 ps enables many modes to be translated with nearly similar efficiencies, as 
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Fig. 11. Parameter study relating to frequency translation in the case with long 
pumps (1000 ps) having peak powers of 400 mW and fiber length 20 m. (a) The 
maximal translation efficiency achieved in the fiber as a function of signal pulse 
width, (b) The length at which the maximal translation efficiency is achieved as 
a function of signal pulse width. 
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Fig. 12. Parametric study relating to the case with short pumps (70 ps) having 
peak powers of 400 mW and fiber length 20 m. (a) The maximal translation 
efficiency achieved in the fiber as a function of signal pulse width, (b) The length 
at which the maximal translation efficiency is achieved as a function of signal 
pulse width. 



can be seen from[T5tc). For this reason we describe this case as "non-discriminatory" 
with respect to its two-color HOM interference properties, since many different types 
of modes lead to good HOM interference. 

Figure ll~4l shows the HOM singular values, as denned by (fl7|) . as a function of Schmidt 
mode number, and the value of P\i for the first three input Schmidt modes as a function 
of fiber length. We emphasize that a value Pi\=0 means that two green or two blue 
photons exit the fiber, but never one green and one blue. This is the HOM interference 
effect as it applies to interference of photons of different colors [10]. Figure [Ml a) shows 
clearly why this case is non-discriminatory in regards to its HOM interference properties; 
many orders of Schmidt modes exhibit HOM singular values near one. 

Figure [T51 shows the absolute amplitudes and phase in both time and frequency of the 
first three green input Schmidt modes (relating to the V matrix). The input amplitudes 
of the pump fields are also shown as a black dashed line. Since the V matrix is nearly 
diagonal, the Schmidt mode amplitudes are almost quantitatively identical to the HG 
basis functions, which stands in contrast to the case in section f4 . 1 1 which had the same 
duration pulse length but twice the pump power. Similarly to the previous cases, the 
phase in time is parabolic with positive curvature, showing that these inputs are also 
down-chirped. 

4-4- HOM interference, short pulse 

Figures [I6j[l7l and [18] detail the case where the pump pulses were 70 ps in duration and 
had peak powers of 200 mW, which is close to the optimum for good HOM interference. 
The goal in studying this case is to contrast it with the former case that had pumps 
with the same peak power but were much longer in duration. This will highlight the 
relative importance of having spectrally broad and temporally short pump pulses in 
regards to the HOM interference behavior. It was found that the optimal characteristic 
time of the HG basis for this case was sa 38 ps. Figure fTBT a) shows the absolute value 
of the V matrix from the SVD of the G g b Green function, while Fig. H7JT b) shows the 
absolute value of the w matrix. The horizontal and vertical axes are the same as they 
were in the previous cases. 

Like the previous cases, both the w and V matrices are nearly diagonal for low Schmidt 
mode number, but unlike the previous case, they quickly become seemingly random 
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Fig. 13. Singular values, V matrix, and w matrix of G g t, Green function relating to 
the case with long pumps (1000 ps) having peak powers of 200 mW optimized for 
good HOM interference, (a) Absolute value of the V matrix, (b) Absolute value 
of the w matrix, (c) Squared singular values p (red) and X" (blue) index by the 
corresponding Schmidt mode number. 



combinations of HG functions for higher-order Schmidt modes. The reason for this is 
the same as it was in the other case having short pump pulses; there are many Schmidt 
modes that have low singular values, and practically any combination of higher-order 
HG functions will have a low, but non-zero singular value. As shown in Fig. ll6f c). only 
the first few Schmidt modes have a singular value meaningfully above zero and will 
therefore have a more ordered mode structure. 

Figure [171 plots the HOM singular values, as defined by (|17|) . as a function of Schmidt 
mode and the value of the two-photon coincidence count probability P\ \ for the first three 
input Schmidt modes as a function of fiber length. Figure [17|a) shows that, in contrast 
to the previous case, the HOM singular values decrease rapidly with Schmidt mode 
number, making this case highly "discriminatory" in regards to its HOM interference 



(a) ' 




(b) ' 



Schmidt mode 




Length [m] 



Fig. 14. HOM singular values 2p n z n and P\\ values relating to the case with long 
pumps (1000 ps) having peak powers of 200 mW. (a) HOM singular value index 
by Schmidt mode number, (b) Two-photon coincidence count probability P\\ as 
a function of length for the first three input Schmidt modes. The red, green, and 
blue lines relate to the first, second, and third Schmidt modes, respectively. 
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Fig. 15. Amplitude and phase in both frequency and time of the first three green 
673-nm input Schmidt modes relating to the case with long pumps (1000 ps) 
having peak powers of 200 mW. The red, green, and blue lines relate to the first, 
second, and third Schmidt modes, respectively. 



properties; only a few Schmidt modes exhibit HOM singular values significantly above 
zero. 



Figure [TBI shows the absolute amplitudes and phase in both time and frequency of the 
first three green input Schmidt modes (relating to the V matrix). The input amplitudes 
of the pump fields are also shown as a black dashed line. Similarly to the previous 
case, these first few Schmidt modes are qualitatively similar to the HG basis functions, 
although they exhibit more asymmetry relative to their width than do the Schmidt 
modes in the previous case. As can be seen in Fig. [18] (c), the frequency amplitudes 
have widths that are comparable to (and, for the second and third mode, somewhat 
larger than) the translation bandwidth, which contributes to the fact that the higher- 
order modes experience a low translation efficiency The phase plots are similar to what 
they were in the previous cases; the phases in time are mostly parabolic with positive 
curvature, implying the inputs are down-chirped. 

Figure [19] shows a parameter study where the pump pulses arc identical Gaussians 
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Fig. 16. Singular values, V matrix, and w matrix of G g t, Green function relating to 
the case with short pumps (70 ps) having peak powers of 200 mW, optimized for 
good HOM interference, (a) Absolute value of the V matrix, (b) Absolute value 
of the w matrix, (c) Squared singular values p (red) and X (blue) indexed by 
the corresponding Schmidt mode number. 
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Fig. 17. HOM singular values and P\\ relating to the case with short pumps 
(70 ps) having peak powers of 200 mW. (a) HOM singular values indexed by 
Schmidt mode number, (b) Two-photon coincidence count probability P\\ values 
as a function of length along the fiber for the first three input Schmidt modes. The 
red, green, and blue lines relate to the first, second, and third Schmidt modes, 
respectively. 
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Fig. 18. Amplitude and phase in both frequency and time of the green 673-nm 
input Schmidt modes relating to the case with short pumps (70 ps) having peak 
powers of 200 mW. The red, green, and blue lines relate to the first, second, and 
third Schmidt modes, respectively. 
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Fig. 19. Parametric study relating to the case with long pumps (1000 ps) having 
peak powers of 200 mW. (a) The minimal P\\ value achieved in the fiber as a 
function of signal pulse duration, (b) The length at which the minimal P\\ value 
is achieved as a function of signal pulse duration. 
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Fig. 20. Parametric study relating to the case with short pumps (70 ps) having 
peak powers of 200 mW. (a) The minimal P\\ value achieved in the fiber as a 
function of signal pulse duration, (b) The length at which the minimal P\ \ value 
is achieved as a function of signal pulse duration. 



with temporal FWHM durations of 1000 ps and peak powers of 200 mW. In Fig. QT] 
(a), the minimal P\\ value achieved along the fiber is plotted as a function of input 
signal pulse FWHM duration. For short signal durations the minimum Pi i value hardly 
changes from 1, which is the initial value (z = 0) of P\\ for all modes, but for long signal 
durations the minimum P\\ value drops to nearly zero. This is most likely due to the 
short signals having large spectral bandwidths and the long signals having small spectral 
bandwidths as compared to the translation bandwidth. This follows because pulses with 
large relative bandwidth experience poor translation as compared to pulses with small 
relative bandwidths. In Fig. [TTTb'). the length at which the minimal P\\ value occurs 
is plotted as a function of input signal duration. For short signal widths, the optimal 
length is the same as the fiber length. This is likely due to the difficulty translating 
signals with large bandwidths; it takes the entire length of the fiber to translate the 
small amount of the signal that is within the translation phase-matching bandwidth. 
This dynamic changes once the signal width is on the order of the phase-matching 
bandwidth, at which point the BS process is no longer limited by bandwidth but by 
the shape and power of the pumps, hence the longer optimal lengths for longer signal 
pulses. 



Figure [20] shows a parameter study where the pump pulses arc identical Gaussians 
with temporal FWHM durations of 70 ps and peak powers of 200 mW. The value of 
J°il, shown in part (a), for short signal inputs behaves similarly to the previous long- 
pumps case: a sharp increase in P\\ for very short signal widths due to the signal having 
larger translation phase-matching bandwidth than the fiber. Unlike the previous case, 
Pn reaches a minimum around 40 ps and begins to slowly increase for wider input 
signal widths. The increase for wider signal widths is likely due to signals beginning 
or evolving to be wider than the pumps, so that only part of the signals wavepackct 
could be translated. The length at which the minimum P\\ occurs is shown in part 
(b). This is a qualitatively similar result to the previous long-pumps cases. For short 
signals the optimal length is equal to the length of fiber (fixed at 20 m), due to the 
excessive bandwidth of the signal in relation to the phase-matching bandwidth. But 
when the signal durations are of the same order as the pumps, the optimal length is 
not as constrained by the phase-matching bandwidth in the same manner, and becomes 
somewhat shorter than the length of the fiber. But when the signals become as long as or 
longer than the pump, the overlap between all the fields in time significantly decreases, 
leading to less translation and raising the value of P\ i . In this regime the fiber length 
also limits the translation, hence the reason the optimal length is equal to the length 
of the fiber. This suggests that were the fiber longer the value of P\\ would be lower, 
which indeed is the case. 

5. Analytic derivation of Schmidt modes 

In order to gain further insight into the numerical simulations presented above, we devel- 
oped an analytical pcrturbativc solution, valid for not-too-high conversion efficiencies. 
The main insights provided by this solution are: (1) The Schmidt modes are found, to 
good approximation, to equal HG functions in this regime. (2) The origin of the input 
Schmidt mode temporal delays seen in Figs. [15] and [18] can be understood as arising 
from differences of group velocities and the requirement for maximal interaction with 
the pump pulses in the fiber. 

The detailed derivation is given in Appendix B, the results of which are summarized 
here. Equations (fT5|) — (fT5|) show how, in general, to express the backward Green functions 
as sums of products of the input Schmidt modes V„ and W n , the output modes v n and w„, 
and the Schmidt coefficients p„ and T„ = (1 — p 2 ) 1 ' 2 . Furthermore, the diagonal Green 
functions G gg and Gtb can be deduced from the off-diagonal Green functions G g b and 
Gbg (and vice versa). In Appendix B we rewrite the third and fourth lines of (|25j) in a 
frame moving with the average group speed of the signals. We include the effects of signal 
convection and second-order dispersion, but do not include the effects of time-dependent 
cross-phase modulation, which were discussed qualitatively after Fig. [7] In this frame 
the dispersion relations for the green and blue signal fields k g (co g ) = fii(O g + /32C0 2 /2 and 
kb((Db) = —fixWb + fiiWb/l- give the associated wavenumbers, correct to second order in 
frequency difference. In our simulations j3j is negative, so the green mode is faster than 
the blue mode. For the parameters in our simulations, the temporal walk-off of pulses 
is significant, but pulse broadening by linear group-velocity dispersion is negligible. 
Therefore, we set y3o = here for the simplest analytical solution. We assume the pumps 
have identical temporal intensity profiles, proportional to exp(— f 2 /2o" 2 ), where a is a 
width parameter. We then find that the perturbative solution for the signals (to first 
order in the coupling parameter 7) involves Green functions that can be expressed in 
SVD form, as are those in (fT3|) (flU)) . In this perturbation theory the Schmidt modes 
are independent of the coupling parameter 7 and pump peak intensities, but of course 



this will break down at higher values of y or pump peak intensities. Notably, we find 
that the Schmidt modes are given by HG functions in this regime. The explicit forms 
of the input green and blue Schmidt modes are, respectively: 



V((o g ) 
W(co b ) 



M 1 


,1/2 



<//„OoC0 ? )exp( 
y/„(foCO fo )exp( 



-icOgfrL/2) 

-ico b lliL/2) 



(32) 
(33) 



where the y/„(foG%) arc the HG functions defined in (26), and the temporal scale fac- 
tor is ?o = (0.621 ct/JiL) ' . Likewise, the output green and blue Schmidt modes are, 
respectively: 



v(a>g) = tj i^(f %)exp(-/© ? j3iL/2), 



w(co b ) 



.1/2 



y/„(t (O b )exp(+i(o b PiL/2). 



(34) 
(35) 



Note that the spectral phases in these mode solutions correspond to temporal delays or 
advances, which in this simple model result in both green and blue modes maximally 
overlapping with the pump pulses at the midpoint (z = L/2) of the medium. In our 
simulations j3i is negative, so the green input mode is delayed and the blue input mode 
is advanced (although note that in our simulations the two pumps walk off from each 
other, which is not the case for the analytic solution given here). 

Using the same parameters as in the numerical simulations leading to Figs. [15] and 
I18[ the approximate input Schmidt modes for the blue signal are plotted as functions 
of time and frequency in Figs. [5T] and [22J (Note that the time-domain modes are also 
described by HG functions.) For conversion efficiencies up to 50%, the approximate 
analytical solutions predict the modes' widths accurately, and their delays or advances 
qualitatively. For higher conversion efficiencies, accurate analytical solutions are not 
yet known. Clearly, the perturbative solutions will not hold for conversion efficiencies 
approaching 100%, as we see significant alterations of the modes in Figs. [71 and [TU1 
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Fig. 21. Approximate Schmidt modes for the faster signal plotted as functions of 
time and frequency, for long pump pulses. The red, green, and blue lines denote 
the first, second and third input Schmidt modes, respectively. 



6. Discussion and conclusions 

In this paper we studied quantum frequency translation (QFT) by the four- wave mixing 
interaction of Bragg scattering (BS) using a realistic numerical model of photonic crystal 
fiber and an analytical approximate solution valid for low conversion efficiency Because 
BS implements a quantum operator transformation analogous to that of an ordinary 
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Fig. 22. Approximate Schmidt modes for the faster signal plotted as functions of 
time and frequency, for short pump pulses. The red, green, and blue lines denote 
the first, second and third input Schmidt modes, respectively. 



beam splitter (mode coupler), it is noiseless (background free) [3J. Furthermore, the 
beam-splitter nature of BS leads to linear-optical quantum operations between weak 
light pulses of different central frequencies, in particular two-photon Hong-Ou-Mandel 
(HOM) interference [TD] . The main questions addressed were: (1) Under what conditions 
does high-efficiency QFT take place? (2) Under what conditions does high-efficiency 
HOM interference take place? (3) Given the shapes of the pump pulses, how can one 
design ideal signal pulse shapes to optimize the above-mentioned processes? (4) How 
are these ideal signal shapes affected by the level of the conversion efficiency for the 
QFT? That is, are the ideal shapes invariants for a given pump shape, or do they evolve 
differently at higher conversion efficiency than at low conversion efficiency? 

Our modeling included realistic input pump and signal fields and effects such as 
convection, dispersion, self-phase modulation and cross-phase modulation. This work 
extends previous theoretical and numerical work [3l HOj on these subjects not only 
through explicitly modeling of specific cases (specific input signal wavepackets), but 
more generally through the determination of the Green function that governs the BS 
process for an entire range of input signal wave packets. Knowing the Green function 
for a process enables the calculation of the signal input and output Schmidt modes and 
corresponding singular values (generalizations of cigenmodes and eigenvalues, respec- 
tively), found by the singular value decomposition of the Green function. Knowledge 
of the Schmidt modes and their associated singular values enables the calculation of 
the effect of the BS process on the input signals without the need to further solve the 
pulse propagation equations. Most importantly for single-photon inputs, this includes 
the frequency translation and two-color two-photon interference effects. It also allows 
for the optimization of these effects. 

In particular, conditions for achieving high translation efficiency of input signal wave 
packets with Gaussian pumps with both long (1000 ps) and short (70 ps) durations at a 
modest peak power of 400 mW were discovered by finding the Green function, and then 
Schmidt modes, of the process. These two cases represent the quasi-CW regime and 
pulsed regime, respectively, in regards to the translation phase-matching bandwidth of 
the process. The long-pumps case was found to be "non-discriminatory," where many 
Schmidt modes were translated with high efficiency. Therefore, this configuration of fiber 
and pumps would be useful for effective translation of many different input signal wave 
packets or in the case where an input packet is a superposition of the high-efficiency 
ones. In contrast, the short-pump case was found to be "discriminatory," where only the 
first few Schmidt modes experienced significant translation. 



Configurations of the discriminatory type would be most suited to applications where 
it is desirable for only a few types of wave packets (ideally only one) to experience 
significant translation, such as in a temporal- mode-selective add/drop filter [16j Hi]- 
Recently it was proposed that such a filter, or "pulse gate" can be implemented using 
frequency upconversion by sum-frequency generation in an optical crystal driven by a 
single strong pump pulse [15]. This process is analogous to the BS process studied here, 
with a major difference being that in upconversion the two signal fields being connected 
by frequency translation must be separated by an amount equal to the pump frequency. 
This limits the process to large frequency separations. The numerical scheme developed 
here is well suited for modeling such a process. 

To gain insight into the FWM Bragg scattering process, we developed an analytical 
perturbation solution. We find that for low efficiency QFT the Schmidt modes are well 
approximated by a set of orthogonal Hcrmite-Gaussian functions with characteristic 
temporal widths determined simply by group- velocity differences and the spectral width 
of the pumps. Nevertheless, we emphasize that the HG modes become strongly altered 
for high conversion, in apparent contradiction to what is argued in |15j . 

Also studied were the conditions for achieving a highly visible two-color photon in- 
terference. We studied two pump configurations having long (1000 ps) and short (70 
ps) durations and a peak power (200 mW) close to the optimum value for a 20-m-long 
fiber. Similar to the high-efficiency translation cases, the long-pump case was "non- 
discriminatory" and the short-pump cases was "discriminatory" in regards to high vis- 
ibility of the two-color photon interference effect. We verified numerically that for the 
long-pumps configuration up to about 50 percent conversion efficiency, both the input 
and output Schmidt mode sets are very well approximated by orthogonal Hermite- 
Gauss (HG) function sets. For the short-pumps configuration, the HG functions are 
reasonable approximations. The approximate theory on which the HG functions are 
based allows one to predict the Schmidt-mode widths and singular values as functions 
of the physical parameters. 

In conclusion, QFT by FWM Bragg scattering appears to be a powerful and robust 
technique for linearly interacting two quantum fields of distinct frequencies. Besides 
offering the ability to change the colors of single photons, it may offer other useful 
abilities in quantum information science, such as linear-optical quantum computing 
"over the rainbow" of optical and near-IR frequencies. 
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Appendix A: probabilities P20 and P02 

The probability for the occurrence of the |2,0) and |0,2) states requires care to derive 
due to the possibility that degenerate photons would be created in the green and blue 
frequency regions. In this case the primed variables attached to the operators cannot 
be taken over to their unprimed counterparts (as co b to a>b was for example). This 
ambiguity can be avoided by calculating P20 directly as the product of the appropriate 
bra and ket vectors (similarly for P02). From ([20]) the term corresponding to the creation 



(36) 



of two green photons leads to an expression for P20 of 

^20 = (V20IIV20) = (vac\ (J J da) g dG)' g a g {G)g)ag{a)' g )A* gg (L,a) g )A* b (L,a) g ) 
x (I I d(O g d(O g a\((O g )a\{(O g )A gg {L,(O g )A gb (L,(O g )\ \vac) 

= (vac\ J J J J dagda'gdagdco'gAlgiL, a> g )A* gb m a)' g )A gg (L, co g )A gb {L, co' g ) 

a g ((Og)a g ((o' g )al(o) g )al((o' g )\vac), 

where the variables of the conjugate wavcfunction are denoted with an overhead bar. 
The key to evaluating this expression is to utilize the field commutation relation. In the 
frequency variable the commutation relation is [a((o),a^((0 )] = 8((Q — (0 ). Hence the 
commutation relation to be used on the innermost operators is 

ag(<B g )a*(o%) = S(cb g - (o g )+a T g ((o g )a g (a) g ). (37) 

This breaks the expression into two parts, and leads to the expression, employing dummy 
variable notation, 



P20 



dcoA gg (L,co)A gb (L,eo) 



dcodco \A gg {L,eo)\ 2 \A gb {L,co )| 2 . (38) 



The expression for Pq2 is the same when A gg — > A bb and A gb — > A bg . The two terms of 
(|38|) have distinct physical interpretations. The first term, derived in part from the delta 
function of Q37p. is essentially the overlap of the two spectral functions A gg and A gb , 
and accounts for cases in which the two green photons are created in the same mode. 
This purely quantum mechanical term arises from the bosonic nature of photons, the 
operators of which obey the above-mentioned commutation relation. The second term 
accounts for the cases in which the two green fields are purely classical, as if these fields 
did not obey the quantum bosonic commutation relations. 

Appendix B: derivation of the approximate Schmidt modes 

According to ([13")) (|T7j|) . the (backward) constituent Green functions are determined by 
the input Schmidt modes V n and W n , the output modes v n and w„, and the Schmidt 
coefficients p„ and T„ = (1 — p 2 ) 1 ' 2 . Furthermore, the diagonal Green functions G gg and 
G bb can be deduced from the off-diagonal Green functions G gb and G bg (and vice versa). 
In this appendix, the off-diagonal Green functions and the associated Schmidt modes 
will be determined approximately, for low conversion efficiencies. 

The third and fourth of (24) govern the signal evolution in the time domain. In the 
frequency domain, the signal equations are 



d z A g (z,(O g ) = ik g (co g )A g {z,co g ) + ijY(z,co g -co b )A b (z,co b )dco bl (39) 

d z A b {z,co b ) = ik b {a) b )A b {z,co b ) + i Y(z,(o b -co g )A g {z,(a g )da) g , (40) 



where C0j is an envelope frequency, k g (o) g ) = jSiO^ + p^O 2 ^ andk b (co b ) = -p\(O b +l52(D b /2 
are the associated wavenumbers (correct to second order in frequency), and y(z,(o) is 
the Fourier transform of the coupling term 2jkA* (z,t)A tj (z,t) divided by 2k. Equations 



(|39|) and (|40|) are valid in a frame moving with the average group speed of the signals 
(so ±j3i are the relative group slownesses) . They include the effects of convection and 
second-order dispersion, but do not include the effect of time-dependent cross-phase 
modulation (which chirps the signals). 

Let Aj(z,(Oj) = Bj(z, (Oj) exp[ikj(cOj)z\. Then the transformed amplitudes Bj obey the 
transformed equations 



d z B g {z,(o g ) = i / y(z,(%-© 6 )exp[%(ffl fe )z-^(©g)z]B fc (z,fi) fe )dfi) fe , (41) 

d z B b (z,(O b ) = i Y(z,(Ob-o) g )e\p[ik g ((o g )z-ik b ((Ob)z}B g (z,(O g )d(Og. (42) 

In general, the B-equations are no simpler than the A-equations, because they depend 
explicitly on z. However, in the low-conversion-efRciency regime, one can replace the 
mode amplitudes on the right sides by the input amplitudes. In this regime, 

B g {L,co g ) w B g (0,(O g )+i / K(z,(O g ,(O b )B b (Q,(o b )dzd(o b , (43) 

B b (L,(o b ) ps B b (0,co b ) + i / K(z,cOb,(Og)Bg(0,(Og)dzd(O g , (44) 

where k(z, (O g ,(O b ) = y(z,(O g — (O b )exp[ik b ((O b )z — ik g ((O g )z) is the kernel in (J4TJ) . Note that 
f{z,CO) = [y(z,-CO)]*, SO K{z,CO b ,(Og) = K*{z,(O g ,(O b ). 

If the pumps are Gaussian in time and identical, and their relative convection (walk- 
off) is neglected, then A p A q = poexp(—t 2 /T 2 ), where po is the peak pump power and T 
is the (Gaussian) pump width. In this case, 

y(co) = [rp cj(2/7r) 1 / 2 ]exp(-o- 2 0) 2 /2), (45) 

where the width parameter o" = t/2 ' . Define the wavenumber- mismatch function 
8((O g , (O b ) = kg((Og) - k b ((Q b ). Then 

S{(Og, co b ) = ft (co g + co b ) + j3 2 (a) 2 - co 2 b )/2. (46) 

By combining the preceding results, and using the identity f e~ z dz = 
Le~ L ' 2 sinc(5L/2), one finds that the integrated kernel 

K(co g , co b ) = [yp La{2ln) l l 2 ] exp[-/j3iL(q ? + co b )/2 - /j3 2 L(o) 2 - (0 2 b )/A\ 

x exp[-o((Og - (O b ) 2 /2]smc[j5iL(co g + co b )/2 + j3 2 L(co 2 - ffl|)/4]. (47) 

The (Oj terms in the first exponential produce delays (advances) in the time domain, 
whereas the ft) 2 terms produce chirps in the frequency domain (convolutions in the 
time domain). By making the approximation sinc(jc) w exp(— cx 2 /2), where c = 0.3858, 
and omitting dispersion (setting p\ = 0) , which is negligible for the parameters of our 
simulations, one finds that 

K{(o g ,(O b ) w [ypoLo(2/7z) l/2 }exp[-il5iL((o g + (O b )/2} 

x exp[-a 2 ((O g - (O b ) 2 /2 - p 2 {(o g + (O b ) 2 /2], (48) 

where the walk-off parameter j3 = c 1 ' 2 j8iL/2. 



Because kernel (j48j) is a Gaussian function of frequency, it can be decomposed using 
the Mehler identity J33] [34] 



exp 



{l+li 2 ){x 2 +y 2 ) 2jj,xy 



2(1 -M 2 ) 



= {n(l-n 2 )} 1/2 Z\n\ n Yn(x)Yn(y), (49) 



where the Hermite- Gaussian functions y/„ were defined in (27). By comparing (|48|) and 
(|4"5|) . one finds that /x = (a — P)/(a+P),x = t g CO g and y = tb<x>b, where ?g = (2j3a) 1//2 = f/,. 
Hence, kernel (|48|) has the singular value (Schmidt) decomposition 

K(co g ,co h ) = ]T A„^(cOg)^(co fo )exp[-/j8iL(ft) i; + ft) i )/2], (50) 

«=o 

where the Schmidt coefficients X n = y/?oL[(cr/j3)(l — n 2 )] l ' 2 \ji\" and the (normalized) 

1/2 

Schmidt modes § n ((Oj) — h \j/„{tjCOj). The kernel is separable when jj, —0 (a = j3), in 

which case Ao = ypo^ and £/ = 2 1 ' 2 <7 = T. In this case, only one mode will be partially 
translated; the others will be unaffected, creating a mode- selective filter, similar to that 
proposed in [15]. 

Let H g i,((O g ,(Ot,) be the Green function that describes the effect on the output trans- 
formed amplitude B g (L,co g ) of the input transformed amplitude Bb(0,(Ob). Then (|4"3"]l 
implies that H ? i(fi) ? ,fi)fc) = iK^^^COi,)- The associated (forward) Green function for the 
original amplitudes, A g (L,(O s ) and Aj(0, ©&), is J g b((Og,(Ot,) — H g b(co g ,(Bt,)exp(if}iL(Dg). 
Hence, 

7^(0%, tofe) = / £ A„0 H ((o g )0„(o)fc)exp[jj3iL(a) s - fl%)/2], (51) 

«=o 

from which it follows that the output (green) and input (blue) Schmidt mode functions 
are given by 

v(cOg) = 0„(to g )exp(-j/3iLco i? /2), (52) 

W(a)b) = ^(0Jfe)exp(-i73iLc>i/2), (53) 

respectively. Likewise, (|4^]) implies that the transformed Green function Hb g ((Ob,(O g ) = 
iK*(co g ,0)b); and the original Green function J bg ((O b ,a) g ) = Hb g ((Ob,(a g )exp(— ifi\L(Db)- 
Hence, 

CO 

J bg {(Ob,co g ) = /^ A„0„((Ofc)0 n ((Og)exp[-j/3iL((o i - co g )/2], (54) 

«=o 

from which it follows that the output (blue) and input (green) Schmidt mode functions 
are given by 

w(a>b) = 0„(a) fc )exp(+i/3iL<B 6 /2), (55) 

V(co g ) = <j> n (co g )exv(+iPiLco g /2), (56) 

respectively. To obtain these, we needed to account for the distinction between forward 
and backward Green functions, and note that here we are propagating A, which corre- 
sponds to the quantum operator a, not a' . If /3i is positive, the green mode is slower 
than the blue mode. The preceding results show that the green input is advanced and 
the blue input is delayed, so their collision is centered on the midpoint of the fiber. This 
location maximizes the distance over which the signals interact with the peaks of the 
pump pulses. In our simulations j3\ is negative, so the green input is delayed and blue 
input is advanced. 



